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Introduction 


The  acronym  LADAR  is  used  here  to  distinguish  optical  radar  systems  that  are  used  to 
detect  and  range  hard-body  targets  from  LIDAR  systems  that  detect  scattered  laser  light 
from  vaiious  diffuse  molecular  and  particulate  constituents  in  the  atmosphere.  Here  we 
consider  a  direct-detection  (monostatic)  LADAR  system  that  emits  a  burst  of  short  (--10  to 
20  ns  Q-switched)  laser  pulses.  These  pulses  are  assumed  to  reflect  off  a  diffuse  target  and 
are  subsequently  collected  by  the  L.ADAR,  and  a  range  determination  of  the  taiget  is 
obtained.  Since  the  forms  available  in  the  literature  for  the  detection  probability  per  pulse 
are  often  untenable  for  practical  systems  studies,  it  is  of  particular  concern  here  to  obtain 
elementary  analytic  approximations  to  the  probabilities  of  detection  in  various 
circumstances  that  are  suitable  for  use  in  commercially  available  spreadsheet  programs 
(c.g.,  E.XCEL,  JAVELIN,  or  123).  Although  the  primary  purpose  of  this  report  is  to 
obtain  results  for  the  probability  of  detection  per  pulse,  we  also  present  an  example  of  a 
target  detection  algorithm,  constructed  from  multiple  pulse  statistics,  which  can  be  used  to 
increase  the  overall  detection  probability  when  the  probabiiitv  ot  detection  |)cr  pulse.  Pj.  is 
iow 


We  consider  targets  that  can  be  both  in  the  near  and  far  field  of  the  receiving  aperture  and 
LADAR  systems  operating  in  the  small  and  large  photocount  regime.  Here  we  assume  for 
precision  that  the  laser  beam  sjxm  size  at  the  target  is  large  in  comparison  to  the  size  of  the 
target  (i.e..  the  transmitted  beam  does  not  resolve  the  target)  and  the  target  depth  fin  the 
range  direction)  is  small  in  comparison  to  the  (longitudinal)  coherence  length  of  the 
transmitted  beam.  When  light,  having  a  deterministic  variation  over  space  and  time  (e  g., 
specular  reflection  from  a  corner  cube  or  a  glint),  impinges  on  a  (square-law) 
photodetector  it  can  lie  shown  that  the  resulting  photoclectron  counts  obey  Poisson 


1 


statistics.  That  is,  we  assume  here  that  the  probability  of  observing  K  photoevents  over  a 
time  period  that  is  long  compared  to  the  optical  period  is  given  by 


P(K)  =  ^e-^',  K  =  0,l,2....  .  (1) 

Jx! 

where  Kg  is  the  mean  number  of  photoevents  expected  per  pulse.  *  This  quantity  is  given 
by 


K. 


(2) 


where  Eg  is  the  mean  return  integrated  signal  fluence  per  pulse  (i.e.,  collected  return 
energy  per  pulse),  h  is  Planck’s  constant,  v  is  the  optical  frequency,  and  q  is  the  detector 
quantum  efficiency.  An  explicit  expression  for  Eg  is  obtained  from  the  LADAR  range 
equation  for  a  diffuse  Lambertian  target  and  a  LADAR  system  with  no  boresight  error  or 
jitter:  2 


E.  - 


-spot .- 


pA, 


n 


T' 

R-y 


where 

Et  =  Transmitted  laser  energy  per  pulse 
^spot  =  Equivalent  laser  beam  spot  area  at  the  target 
p  =  Diffuse  target  refectance 
Ax  =  Target  area,  normal  to  the  line-of-sight 
R  =  Propagation  range 
Ar  =  Receiver  clcar-aperture  area 
T  =  One-way  propagation  path  transmittance 


(3) 
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For  propagation  through  vacuum,  T  is  just  the  LADAR  system  optics  transmittance.  Here, 
we  assume  a  laser  beam  irradiance  profile  given  by 


I(6^  =  I<,exp[-  <9^/24] 


(4) 


where  Iq  is  the  beam  irradiance  for  0=0.  For  systems  that  employ  a  fixed  beam 
divergence,  we  thus  have 


A,^„.=2;r(R^,)^ 


(5) 


We  note  that  the  l/V7  full-angle  beam  divergence  is  equal  to  20q,  and  the  corresponding 
1/^2  full-angle  beam  divergence  equals  2(20o)  =  400-  Thus,  for  example,  a  beam 
characterized  here  by  Oq  =  5prad  corresponds  to  a  beam  with  a  20|.irad  He-  full -angle 
beam  divergence.  Such  a  system  results  in  energy  returns  that  are  proportional  to  R"^.  On 
the  other  hand,  LADAR  systems  employing  appropriate  zoom  optics  that  keep  = 
constant  at  the  target,  independent  of  range,  result  in  energy  returns  proportional  to  R'-. 
Here  we  assume  uniformly  illuminated  targets  for  simplicity.  The  case  of  nonuniformly 
illuminated  targets  can  be  formally  included  by  replacing  Aj  in  the  LADAR  range  equation 
by  some  effective  target  area  (Aj)^,Y  (<  A-j-).  For  example,  for  a  target  possessing  a 
circular  cross-sectional  area  normal  to  the  line-of-sight  and  a  LADAR  employing  a  laser 
beam  spread  of  Gaussian  shape,  it  can  be  shown  that  (Aj)gff  =  nryf,  where  the  reduction 


factor  f  =  [2  6ll6\) 


I-.J.  is  the  radius  of  the  target,  Oj 


rj/R,  and  0o  is  the 


\|^fe  angular  laser  beam  radius.  Uniform  illumination  is  obtained  for  0j  «  0o,  where 
f  =  1  and  (Ax)eff  “  Aj.  Although  very  few  optical  radar  transmitters  resolve  the  target, 
this  case  can  be  incorporated  into  the  present  analysis  simply  by  setting  Aj  =  A^po,  in 


Eq.  (3). 
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In  many  cases  of  practical  concern,  the  return  signal  is  not  deterministic,  but  rather  has  a 
random  variation  over  space  and  time  (e.g.,  from  pulse  to  pulse).  It  is  now  necessary  to 
treat  the  Poisson  distribution  of  Eq.  (1)  as  a  conditional  probability  distribution,  the 
conditioning  being  based  on  knowledge  of  the  probability  distribution  of  the  return 
integrated  (over  the  receiver  collecting  aperture)  signal  fluence.  In  practice,  however,  it  is 
the  unconditional  probability  distribution  of  photocounts  that  relates  to  required  system 
performance  quantities  such  as  the  probability  of  detection.  The  unconditional  probability 
distribution  is  obtained  by  averaging  the  conditional  Poisson  distribution  of  Eq.  (1)  over 
the  statistics  of  the  return  integrated  signal  fluence.  The  details  of  the  calculation  have 
been  documented  in  Reference  (3)  with  the  result  that 


P(K)  = 


r(K  +  m) 
r(K  +  i)r(m') 


•n-K  , 


-f 


m 


(6a) 


E[K]-K, 


(6b) 


Var[Kl=K, +  .  (6c) 

m 

where  E  [•]  and  Var  [•]  denote  the  expectation  value  and  variance,  respectively,  F  ( )  is  the 
gamma  function,  and  m  is  given  by 


m-l  +  ^^^  (7) 

A'  R’ 

The  first  and  second  term  on  the  right-hand  side  of  Eq.  (6c)  is  the  signal  shot-noise  and 
speckle-noise  contribution  to  the  signal  variaiice,  respectively.  As  discussed  in  Reference 
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(3),  m  can  be  interpreted  as  the  mean  number  of  independent  irradiance  correlation  cells 
(i.e.,  speckles)  contained  within  the  collecting  aperture.  Examination  of  Eq.  (6)  reveals 
that  the  signal  statistics  depend  not  only  on  the  mean-signal  photoelectron  count  per  pulse, 
Kg,  but,  in  addition,  also  on  the  mean  number  of  correlation  cells,  m,  (or  speckles) 
contained  within  the  collecting  aperture.  Since  the  signal-count  fluctuations  increase  with 
decreasing  values  of  m,  it  is  seen  that  LADAR  system  performance  generally  deteriorates 
as  the  number  of  speckles  decreases  with  the  worst  performance  obtained  for  m  =  1.“^  In 
addition,  it  is  easily  seen  from  Eq.  (6)  that  the  Poisson  distribution  is  an  accurate 
approximation  to  the  negative-binomial  when  .n  »  Kg.  Thus,  for  an  optically  rough 
target,  the  assumption  of  Poisson  signal  statistics  is  justified  only  when  the  average  number 
of  signal  photoelectrons  per  spatial  correlation  cell  is  much  less  than  unity.  For  any  given 
physical  situation,  both  Kg  and  m  are  obtained  from  Eqs.  ^2)  and  (7).  and  the  results  are 
used  in  Eq.  (6). 

The  distribution  given  by  Eq.  (6)  is  known  as  the  negative-binomial  distribution.  For  m 
near  unity  (unresolved  targets),  a  single  speckle  is  contained  in  the  collecting  aperture,  the 
collected  energy  obeys  negative-exponential  statistics,  and  the  corresponding  photocounts 
obey  Bose-Einstein  statistics.  On  the  other  hand,  for  m  »  1  (resolved  targets),  many 
speckles  are  contained  within  the  collecting  aperture,  and  the  collected  energy  distribution 
behaves  as  a  delta  function  centered  at  the  mean  value  Eg.  Hence,  in  the  limit  m  ->  x, 

P(K)  must  obey  Poisson  statistics.  These  analytic  features  are  contained  in  the  negative- 
binomial  distribution  given  by  Eq.  (6). 

Figure  1  is  a  plot  of  m  as  a  function  of  range  for  a  wavelength  of  0.532|am  and  a 
collecting-aperture  diameter  of  20  cm.  Examination  of  Figure  1  reveals  that  for  R  >  1.5 
Megameters,  m  '>■  1  for  Aj  5  1  m^.  As  a  result,  the  return-signal  photocount  is  well 


approximated  by  Bose-Einstein  statistics.  On  the  other  hand,  for  R  <  0.5  Mm,  m  »  1  for 
Aj  >  5  m2  and  hence  the  signal  photocount  is  well  approximated  by  Poisson  statistics  if 
Kg  «  m. 

We  now  distinguish  between  the  low  signal  photoelectron  count  regime,  where  Kg  «  100, 
and  the  high-count  regime,  where  Kg  >  100.  In  the  low  photocount  regime,  the  discrete 
character  of  the  signal  counts  is  retained  while  for  the  high  photocount  regime  it  is 
customary  to  characterize  the  collected  signal  by  a  continuous  distribution  of  signal  current 
i,  related  to  the  corresponding  signal  photocount  by 


(8) 


where  q  is  the  electronic  charge,  and  x  is  the  reflected  return  laser  signal  pulse  width. 

As  an  example  illustrating  the  expected  value  of  signal  counts  per  pulse  obtained  under 
various  circumstances,  we  consider  the  following  values  of  LADAR  system  parameters: 

Ej  =  0.5  J/pulse,  a  collector  diameter  of  20  cm,  T-  =  0.5,  X  =  0.532tim,  Aj  =  5  m^,  p 
=  0.4,  and  t]  =0.8.  We  consider  a  system  that  employs  a  fixed  laser  beam  width  as  well 
as  one  that  keeps  a  fixed  spot  at  the  target,  independent  of  range.  The  results  for  Kg  as 
obtained  from  Eq.  (2)  are  plotted  in  Figures  2  and  3  as  a  function  of  range  for  various 
values  of  Gq  and  Agp^^,  respectively.  Examination  of  Figure  2  reveals  that,  for  the 
parameters  assumed  above,  the  low  photoelectron-count  regime  (Kg  «  100)  is  expected  for 
R  >  0.4  -  1.3  Mm,  depending  on  the  value  of  Gq  for  the  fixed-beam-divergence  system. 

On  the  other  hand,  examination  of  Figure  3  reveals  that  for  the  fixed-spot-on-target  case 
the  corresponding  low  photocount  regime  is  expected  for  R  >  1.0-3  Mm. 
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Figure  2.  The  mean  number  of  signal  photocounts  as  a  function  of  range  for  various 
values  of  laser  beam  divergence.  Note  that  the  beam  angles  indicated  here 
correspond  to  the  l/Ve  angular  beam  radius.  The  l/tr‘  full  angle  beam 
divergence  is  equal  to  four  times  the  angles  indicated  above. 
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Figure  3.  The  mean  number  of  signal  photocounts  as  a  function  of  range  for  various  values 
of  the  laser  spot  area  at  the  target. 


In  ihe  second  section,  we  develcp  elementary  analytic  expressions,  suitable  for 
commercially  aN'ailabfe  spreadsheet  programs,  such  as  JAVELIN,  for  the  probability  of 
detection  including  the  effects  of  both  unintentional  beam  boresight  error  and  beam  jitter  in 
the  low  pho'oconnt  regime.  In  the  third  section,  we  consider  the  high  photocount  regime 
and  develop  analytical  expressions  for  the  corresponding  probability  of  detection.  Finally, 
in  the  Summary,  we  summarize  the  approximations  of  the  detection  probabilities  obtained 
here,  together  with  their  respective  regions  of  validity.  Those  readers  who  want  to  use 
directly  the  respective  approximations,  but  do  not  want  to  get  involved  with  the  technical 
details,  can  skip  to  the  final  results  in  the  Summary. 

In  order  for  this  report  to  be  as  self-contained  as  possible,  we  have  included  a  number  of 
appendices  to  serve  as  technical  references  and  document  some  interesting  results  that  may 
not  be  well  known.  In  Appendix  A,  we  assume  a  LADAR  system  with  circularly 
symmetric  Gaussian  jitter  with  a  fixed  pointing  offset  and  derive  the  corresponding 
probability  distribution  of  intensity.  In  Appendix  B,  we  show'  how  the  general 
photoelectron  count  distribution,  given  by  the  negative-binomial  probability  distribution, 
reduces  to  the  Gaussian  (i.e.,  normal)  distribution  for  both  large  photocounts  and  speckle 
numbers  contained  within  the  collecting  aperture.  In  addition,  we  show  how  the  negative- 
binomial  distribution  reduces  to  the  Bose-Einstein  and  Poisson  distribution  in  the  limit  that 
the  speckle  number  tends  to  one  and  infinity,  respectively.  In  Appendix  C,  we  consider 
signal-noise-limited  detection  for  large  photoelectron  counts  and  mean  number  of  speckles 
near  unity  and  derive  the  corresponding  count  distribution  and  probability  of  detection.  In 
this  regard,  we  have  derived  a  possible  new  and  interesting  mathematical  result  involving 
the  incomplete  gamma  function.  Finally,  in  Appendix  D,  we  consider  the  case  of  a  signal 
obeying  Bose-Einstein  photocount  statistics  (such  as  that  arising  from  unresolved  targets), 
which  is  immersed  in  Poisson  noise,  and  derive  the  photocount  probability  of  signal  plus 
noise  as  well  as  the  corresponding  probability  of  detection.  In  particular,  we  obtain  a 
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closed-form,  anaJytical  expression  for  the  probability  distribution,  valid  for  all  values  of 
signal  and  noise  photocounts. 
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Probability  of  Detection  for  Low  Photoelectron  Counts 


Detection  Statistics  without  Pointing  Errors 

First  consider  the  sifu'^tion  where  the  mean  photocount  number  is  small,  say  less  than 
about  50.  For  the  example  considered  in  the  previous  section,  we  have  seen  that  the  low 
photocount  regime  is  expected  to  be  obtained  for  ranges  >  1  Mm.  For  a  LADAR  system 
to  be  viable  in  this  regime,  it  is  necessary  to  assume  that  over  a  pulse-length,  the 
background,  and/or  photodetector,  noise  is  much  smaller  than  the  corresp)onding  return 
signal  strength.  That  is.  in  the  low  photocount  regime,  we  assume  a  (near)  quantum-noise- 
iimited  LADAR  system.  Such  physical  circumstances  may  be  obtained  for  a  space-based 
LADAR  system  with  sufficiently  narrow  field-of-vie\v  and  optical  bandpass  that  operates  in 
the  natural  environment  against  a  deep  space  background  (i.e..  no  earth,  moon.  sun.  or 
stars  in  the  field-of-view)  and  with  the  line-of-sight  between  the  LADAR  and  the  target 
having  tangent  altitudes  greater  than  '-100  km.  In  addition,  it  is  necessary  to  employ  a 
sufficiently  low  noise  pholodetector  (e  g.,  a  cooled  photomultiplier  or  a  semiconductor 
avalanche  photodiode  operating  above  breakdown  in  the  triggered  avalanche  Geiger  mode). 

Since  in  the  low  photocount  regime  we  assume  negligible  background  and/or  detector- 
induced  noise  photocounts  over  the  measurement  time  of  interest,  the  probability  of  a 
detection,  P^,  of  a  single  return  signal  pulse  is  equal  to  the  probability  that  at  least  one 
photocount  occurs  during  this  lime  interval.  This  quantity  is  just  equal  to  one  minus  the 
probability  of  obtaining  zero  counts  during  this  time  interval.  Hence,  from  Eq.  (6)  we 
obtain 
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p.= 

where  Kj  is  the  mean  number  of  signal  photoelectron  counts  expected  from  the  LADAR 
pulse,  and  m  is  the  mean  number  of  speckles  of  the  diffuse  reflected  pulse  contained  within 
the  collecting  aperture. 

In  LADAR  systems,  where  the  probability  of  detection  per  pulse,  Pj,  is  low,  target 
detection  algorithms,  constructed  from  multiple  pulse  statistics,  can  be  employed  to 
increase  the  detection  probability.  For  example,  when  N  independent  pulse  measurements 
are  combined,  the  probabiliity  of  detecting  at  least  one  photocount  from  a  burst  of  N  pulses 
is  obtained  as 


1  -  P{0) 


P,(N)=1-(1-P,y'  =  1-P(0)" 


-  1 


Nm 


(10) 


Thus,  by  combining  N  independent  pulses,  the  mean  photocount  Kj  per  pulse  needed  for  a 
specific  detection  probability  is  lowered,  and,  consequently,  for  a  given  Kg,  the  acquisition 
range  is  correspondingly  extended.  Of  course,  this  improvement  is  at  the  cost  of  the 
energy  to  produce  the  added  pulses  and  the  time  required  to  make  the  measurements. 
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Detection  Statistics  with  Pointing  Errors 


In  general,  the  effects  of  pointing  errors  result  in  a  reduction  of  target  illumination  and  a 
corresponding  decrease  in  return-signal  energy.  As  a  result,  the  LADAR  detection 
probability  is  reduced  in  comparison  to  the  corresponding  value  in  the  absence  of  fKtinting 
errors.  In  the  following,  a  fixed  angular  pointing  offset  bias  and  circularly  .symmetric 
Gaussian  jitter  are  assumed. 

The  quantity  P^,  given  by  Eq.  (9),  represents  the  detection  probability  of  a  single  pulse 
without  any  pointing  errors.  The  corresponding  probability  with  beam  jitter  and  offset  is 
obtained  by  integrating  P^j  over  the  mean  photocount  distribution  as  a  function  of  the 
pointing  statistics.  That  is,  with  jitter  and  offset,  the  mean  photocount  is  now  considered 
as  a  random  variable,  k,  with  pointing  statistics  given  by  (see  Appendix  A) 


f'p.(k)  = 


a^exi 


p(-  a'bV'i) 


(VK.f''j„(a^b727n^/K.)).  0<k<K, 


=  0,  otherwise 


(11a) 


where  circularly  symmetric  Gaussian  jitter  statistics  are  assumed,  Jq  is  the  Bessel  function 
of  the  first  kind  of  order  zero,  0p,  is  the  (angular)  jxiinting  offset,  oj  is  the  single-axis 
Gaussian  jitter  standard  deviation,  and  ©q  is  the  Gaussian  laser  beam  shape  standard 
deviation. 


a 


(lib) 
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(11c) 


9p,  is  the  (angular)  pointing  offset,  Oj  is  the  single-axis  Gaussian  jitter  standard  deviation, 
and  @Q  is  the  Gaussian  laser  beam  shape  standard  deviation. 

The  probability  of  detection  of  a  pulse  with  pointing  errors,  denoted  by  Pd ,  is  thus 

K., 

P,  =  j  PAk)PJk)dk  .  (12) 

where  P,j(k)  is  given  by  Eq.  (9),  with  replaced  by  k.  and  Pp|(k)  is  given  by  Eq.  (il). 

In  general,  closed- form  solution  exists  for  this  integral,  although  numerical  results  can  be 
readily  obtained  on  a  digital  computer. 

One  would  not  expect  to  design  and  employ  a  LADAR  system  with  excessive  pointing 
offset  and  jitter.  To  this  end,  we  have  derived  an  elementary  engineering  approximation 
for  Pd  valid  for  a,  /Gq  <  0.5  -  0.75  (or  a  =  Go/a,  >  1.5  •  2).  and  Gp,/0o  <  1.0.  The  most 
accurate  approximation  obtained  so  far  that  is  suitable  for  spreadsheet  programs  is  given  by 


P~i^PAf(^) 


=  1- 


1  + 


K, 


m 


(13) 
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-Ji 


where 


K,=j  lcP,.(k)dk 

0 


(14) 


In  obtaining  the  engineering  expression  given  by  Eq.  (14),  we  expanded  the  Bessel 
function  in  Eq.  (11)  to  terms  of  second  order  in  b,  performed  the  integration  over  k,  and 
then  exponentialized  the  result  so  that  the  limit  of  Ks  as  Oj  -v  0  is  identical  to  the  exact 
answer;  Ks  =  Kg  exp['b2/2]. 


A  comparison  of  the  exact  numencaJ  and  approximate  results  are  shown  in  Figures  4 
through  9  for  a  wide  range  of  the  parameter  values.  We  note  that  any  practical  LADAR 
system  design  would  not  allow  the  laser  beam  width,  0q,  to  be  less  than  both  the  jitter 
standard  deviation  Oj  and  the  offset  bias  8p,.  A  reasonable  design  goal  would  limit  both 
the  ratio  Oj/Go  and  0p,  to  be  less  than  about  1/2  (i.e.,  a  >  2,  and  b  <  0.5).  Examination  of 
Figures  4  through  9  reveals  that  over  a  wide  range  of  the  parameters  m  and  K^,  the 
approximate  analytical  result  of  Eq.  (13)  agrees  within  10  percent  of  the  numerical  result 
for  both  a  >  2  and  b  <  1 .  That  is,  in  the  range  where  a  LADAR  system  should  be 
designed,  Eq.  (13)  represents  a  good  engineering  approximation  for  the  probability  of 
detection  in  the  presence  of  jitter  in  the  low  photocount  regime. 
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b  =  0,  m  =  1.5 


curves  are  the  exact  numerical  results,  and  solid  curves  are  the  approximate 
numerical  results  (given  by  Eq.  13),  respectively. 


a 


Figure  6.  Probability  of  detection  for  the  case  where  m  =  1 .5  and  b  =  1 .  The  dashed  curve 
is  the  numerical  result,  and  solid  curves  are  the  numerical  and  approximate  result 
(given  by  Eq.  13),  respectively. 
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» 


Figure  8.  Probability  of  detection  for  the  case  where  m  =  15  and  b  =  0.5.  The  dashed  curve 
is  the  numerical  result,  and  solid  curves  are  the  numerical  and  approximate  result 
(given  by  Eq.  13),  respectively. 


Probability  of  Detection  for  High  Photoelectron  Counts 


Detection  Statistics  without  Pointing  Errors 

Consider  situations  where  a  large  number  of  both  signal  and  noise  photocounts  per  pulse 
are  expected,  such  as  a  space-based  LADAR  system  that  views  below  the  horizon. 
Following  G-jodman'^  and  Kingston^  we  assume  that  the  total  non-signal-related  noise  (i.e., 
background,  dark  current,  thermal,  etc)  is  Poisson  distributed.  For  large  photocounts,  this 
distribution  is  to  a  very  good  approximation  Gaussian.  As  discussed  previously,  it  is 
customary  in  the  high  photocount  regime  to  deal  with  currents  rather  than  count  number 
[see  Eq.  (8)].  Thus  the  ac  component  of  the  non-signal-related  noise  current  from  the 
detection  system  can  be  expressed  as 


Pn(i)  = 


1  r  i' ' 

-===exp! - r 

yl^TTCx;  L  j 


(15) 


where  the  non-signal  noise  current  variance  is  given  by 


T-  ' 


(16) 


where  Kg  is  the  mean  background  photoelectron  count  per  measurement  interval  T  (t  >  t)  . 
o jj.  and  are  the  dark  current  and  thermal  noise  current  variance,  respectively. 


The  probability  of  false-alarm,  Pfg,  is  obtained  by  integrating  Eq.  (15)  above  a 
measurement  threshold  current  ij.  We  thus  obtain 
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(17) 


Pf,= 


1  -  erf 


(TNR),^ 

■J2  J 


where  erf  denotes  the  error  function,  and  the  threshold  to  (non-signal  related)  noise  ratio  is 
given  by 


(rNR).  (18) 

O'. 

A  curve  fit  to  the  numerical  inversion  of  Eq.  (17).  which  is  accurate  to  better  than  1 
percent  o\'er  the  range  10"  '  ^  Pfa  ^  i0'“,  is  given  by 

(TNR)o  =  0.513  -  1.02  log  Pf,  -  0.0539(log  Pfe)2  .  (19) 

For  example,  a  required  false-alarm  probability  of  10'^  requires  that  the  (TNR)^  3. 1. 

We  note  that  Pfg  is  determined  by  the  non-signal-related  noise. 

The  probability  distribution  of  signal  photocounts  is  given  by  the  negative-binomial 
distribution  of  Eq.  (6).  Now  for  both  »  1  and  m  »  1,  it  can  be  shown  that 
probability  distribution  of  the  signal  current  distribution  is  also  Gaussian  and  given  by  (see 
Appendix  B) 


P.O)  s 


yj23T(Tl 


-exp 


2a;  j 


(20) 
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where  the  mean  signal  current 


(21a) 


and 


(21b) 


The  Gaussian  probability  distribution  is  a  good  approximation  to  the  photoelectron 
distribution  as  long  as  the  standard  deviation  is  small  compared  to  the  mean.  In  practice, 
this  is  well  statisfied  for  Kg  >  50-100  and  m  >  5-10. 

Since  both  the  signal  and  noise  are  Gaussian  variates,  the  probability  distribution  of  both 
signal  plus  noise  is  also  Gaussian  and  given  by 


P,n(i)- 


1 


(22) 


where 


o:=o:  +  <r 


(23) 


The  probability  of  detection  is  obtained  by  integrating  above  a  measurement  threshold 
current  if.  The  result  is  given  by 
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(24) 


Pa  = 


\ 

2 


“ 

Is  -‘t 

■ 

1  -t-  erf 

We  note  that  the  probability  of  detection  is  a  function  of  both  the  signal-to-noise  ratio 
(SNR)  and  the  threshold-to-noise  ratio  (TNR),  where  the  noise  includes  the  contribution 
due  to  signal  fluctuations,  as  is  evident  in  Eq.  (23). 

Next  we  consider  the  case  of  Kg  »  1,  but  m  near  unity  (<  10).  In  this  case,  it  is  more 
useful  to  return  to  the  discrete  counting  distribution.  Goodmar.3  has  shown  that  the 
probability  distribution  of  the  total  number  (signal  plus  noise)  of  phoioelectrons  ejected 
during  reception  of  a  return  pulse  is  given  by 


P,n(K)  = 


m 

rexp[-n],^  r(K-j4-m) 

f  K.  ^ 

Im  -hKJ 

Hm)  r(j  +  i)r(K  -  j  + 1) 

Im  +kJ 

(25) 


where  n  is  the  mean  number  of  non-signal-related  noise  photoelectrons  emitted  during  the 
measurement  period  T  (i.e.,  n  =  T’  a’/q').  In  general,  it  can  be  shown  by  methods 
discussed  in  Chapter  9  of  Reference  3  that  E[K]  =  n  -I-  Kg  and  Var  [K]  =--  n  •+•  Kg  -h 
Klim. 

Although  suitable  for  numerical  computations  on  a  digital  computer,  Eq.  (25)  is  not  useful 
for  spreadsheet  programs  such  as  JAVELIN.  Under  certain  circumstances,  simplification 
of  Eq.  (25)  can  be  obtained.  For  Kg  »  1  and  Kg  »  n  (i.e.,  signal-noise-limited 
detection),  it  can  be  shown  that  (see  Appendix  C) 
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(26) 


which  is  valid  for  m  >  1.  This  is  known  as  the  gamma  distribution.  Converting  to  current 
and  integrating  Eq.  (26)  above  a  measurement  threshold  yields  the  probability  of  detection 


Pa  = 


r(m,g^) 
r(m) 


(27) 


where 


=  m~ 

m 

1, 


(28) 


and  r(a,z)  is  the  complementary  incomplete  gamma  function  given  by^ 


r(a,z)=:  Je-'t“'dt  (29) 


Although  no  analytical  approximation  exists  for  arbitrary  m,  closed-form  expressions  can 
be  obtained  for  integer  m.  Asa  result,  we  present  analytical  results  for  Pj  for  the  first  few 
values  of  integer  m. 

m  =  1; 

=exp[-iT/i.]  ,  (30a) 
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m  =  2; 

P,  +1)  ,  (30b) 

m  =  3; 

Pd=^(s?+2®j+2)  .  (30c) 

m  =  4; 

?,=^(sl+3e;+6s,+6)  ,  (30d) 

m  =  5; 

P,  =  [si  +  4sl  +  I2£l+ 24  e, +24)  ,  (30e) 

m  =  6; 

P,  =  (si  +  5^  +  20£^  +  60£^  ■+120^+120)  ,  (30f) 

For  noninteger  m,  it  is  suggested  to  smoothly  interpolate  Eqs.  (30,  a  through  f)  between 
successive  integers.  For  larger  values  of  m,  the  expressions  based  on  the  Gaussian 
probability  distribution  of  Eq.  (22)  with  =  Og  should  apply  (see  Appendix  C).  Hence 
for  m  >  5-10,  we  can  employ  P^  as  given  by  Eq.  (24). 

Except  for  the  case  m  =  1 ,  we  have  not  been  able  to  obtain  closed-form  expressions  for 
Pgn  for  Kg  »  1,  small  m  and  finite  n.  For  m  =  1,  it  can  be  shown  that  (see  Appendix  D) 
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„  /T.N  expf--(K-n)/K3lr(K  +  l,n) 

K,  r(K  + 1) 

K. 

where  U  is  the  unit  step  function.  The  probability  distribution  of  signal  plus  noise  current 
can  thus  be  expressed  as 


Psn(i)  = 


exp[-(i-ij/i, 


uO-i.) 


(32) 


where  in(=qn/T)  is  the  mean  non-signal-related  noise  current.  The  corresponding 
probability  of  detection  is  now  obtained  as 

P,  ^  e-’  ,  (33) 

where  the  threshold  current  is  now  measured  relative  to  the  noise  floor  current 

Detection  Statistics  with  Pointing  Errors 

Just  as  in  the  low  photoelectron  current  regime,  the  most  accurate  approximation  to  include 
the  effects  of  a  pointing  offset  bias  and  jitter  is  to  replace  the  mean  signal  (ij  or  K^)  by  the 
corresponding  value  obtained  by  averaging  over  the  pointing  error  probability  distribution 
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given  by  Eq.  (11).  That  is,  Kg  (and  4,  since  it  is  proportional  to  the  signal  photoelectron 
count)  is  replaced  by  Ks  (and  the  corresponding  expression  for  is)  given  by  Eq.  (14). 

We  consider  the  regime  where  both  K  and  m  »  1.  In  this  case  Pd ,  the  detection 
probability  averaged  over  the  pointing  statistics,  is  obtained  from  Eq.  (24)  with  ig  replaced 
by  is,  where 


is 


i 


/'  2  \ 

1  «  1 

I 

1 - 

ri 

'  ‘  2  1 

va  +  IJ 

2  ’ 

>^3“  +  1/ 

(34) 


Figures  10  through  12  show  a  comparison  of  the  numerical  and  approximate  result*  for  the 
detection  probabilities  as  a  function  of  a  =  Go/oj  for  a  representative  range  of  parameter 
values.  In  these  figures,  the  (current)  SNR  and  the  corresponding  TNR  are  defined  as 


SNR  = 


(35a) 


and 


TNR  = 


(35b) 


Note,  in  particular,  that  the  noise  variance,  Os„,  includes  signal  fluctuations,  in  contrast  to 
the  false-alarm  probability,  which  is  based  on  non-signal-related  noise  only. 
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Figure  10.  Probability  of  detection  for  a  large  number  of  photocounts  and  m  is  large  (>  5-10) 
and  b  =  0.  The  dashed  and  solid  curves  are  the  numerical  and  approximate  result, 
.respectively. 


Figure  1 1 .  Probability  of  detection  for  a  large  number  of  photocounts  and  m  is  large  5-10) 
and  b  =  0.5.  The  dashed  and  solid  curves  are  the  numerical  and  approximate 
result,  respectively. 


respectively. 


Examination  of  Figures  10  through  12  reveals,  as  in  the  low  photoelectron  count  regime, 
that  the  algorithm  used  here  to  include  the  effects  of  pointing  errors  is  sufficient  for  use  as 
an  engineering  approximation  over  the  design  range  of  a  practical  LADAR  system.  We 
note  that  evaluation  of  the  error  function  may  not  be  compatible  with  some  spreadsheet 
programs.  If  this  is  the  case,  use  of  a  global  rational  approximation  to  the  error  function 
can  be  employed,  such  as®: 


erf(x)  =  I  --  (a,i  +  a,t‘  +  a3t^)e~’‘'  +  ^(x)  , 

where 

jc:U)js2.5  X  lU'-  , 

1 

r  = - , 

1  +  py 


and 


p  =  0.47047,  a,  =  0.34802,  a,  =  -0.09587,  and  a3  =  0.74786  . 

Additionally,  we  consider  the  case  of  signal-limited-noise  detection  where  Kg  »  1,  and  m 
is  near  unity.  In  this  case,  we  have  from  Eqs.  (27)  and  (28) 


Pd 


T{m,Em) 

r(m) 


(36) 
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where 


Sm-m’^=n^(ijis)  ,  (37) 

and  the  ratio  i  $  /  ijis  given  by  Eq.  (34).  For  low  integer  values  of  m,  the  results  given  by 
Eq.  (30)  are  to  be  used  with  £„  replaced  by  Sm- 

Figures  13  through  15  show  a  comparison  of  the  numerical  and  approximate  results  for  the 
detection  probability  as  a  function  of  Go/oj  for  m  =  1  and  a  representative  range  of 
parameters.  The  results  for  m  greater  than  unity  behave  similarly.  Examination  of  these 
figures  reveals  that  the  approximations  obtained  above  are  sufficient  for  engineering 
purposes. 

Finally,  we  note,  that  when  N  multiple  independent  measurements  are  combined,  the 
probability  of  detection  becomes 


P<(N)=I-(1-P,f 


(38) 


where  Pj  is  the  detection  probability  of  a  single  pulse. 
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Figure  13.  Probability  of  detection  for  m  =  1,  b  =  0  and  both  the  SNR  and  TNR  are  large. 
The  dashed  and  solid  curves  are  the  numerical  and  approximate  results, 
respectively. 


Figure  14.  Probability  of  detection  for  m  =  1,  b  =  0.5  and  both  the  SNR  and  TNR  are  large. 
The  dashed  and  solid  curves  are  the  numerical  and  approximate  results, 
respectively. 


Summary 


In  this  report,  we  have  derived  accurate  (to  within  10  percent),  elementary,  analytical 
expressions  for  the  probability  of  detection,  Pd ,  over  a  reasonable  design  goal  range,  of  a 
monostatic,  direct-detection  LADAR  system  in  the  presence  of  unintentional  pointing 
errors.  We  distinguish  between  two  photocount  regimes— iow  (Kg  «  100)  and  high  (Kg 
greater  than  50-100).  As  discussed  in  the  second  section,  in  the  low  photocount  regime, 
we  are  primarily  interested  in  signal-noise-limited  LADARs.  In  this  regime,  we  obtain 
results  for  Pd  that  are  valid  for  resolved  targets  (m  >  10),  partially  resolved  targets  (1  < 
m  <  5-10),  and  unresolved  targets  (m  «  1).  On  the  other  hand,  in  the  high  photocount 
regime,  we  obtain  results  for  Pd ,  valid  for  both  m  »  1  and  m  »  1  and  arbitrary  values  of 
non-signal-related  noise.  For  partially  resolved  targets,  we  have  not  been  able  to  obtain 
elementary  approximations  of  Pd  except  for  cases  of  large  values  of  the  (SNR)^, 

where 


(SNR)„ 


(39) 


is  the  signai-to-(non-signal  related)  noise  ratio. 

In  Table  1,  we  summarize  several  of  the  key  parameters  used  here  with  a  reference  to  their 
appearance  in  the  body  of  the  paper.  For  the  low  photocount  regime  (Kj  «  100),  the 
detection  probability  of  a  burst  of  N  independent  pulses,  averaged  over  the  pointing 
statistics,  is  given  by 
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m>\  , 


(40) 


V  m  ■ 


-mN 


where 


(  1 

■_ir  ba^  ^ 

2" 

2U'  -nj 

(41) 


In  Table  2,  we  summarize  the  results  obtained  for  Pd  in  the  high  photocount  regime 
(Kj  >  50-100).  For  a  burst  of  N  independent  pulses  the  probability  of  detection  is  given  by 
Eq.  (38). 

The  approximations  for  Pd  given  in  Eq.  (40)  and  Table  2  are  accurate  to  within  10  percent 
for  a  =  Oo/Oj  >  2,  and  b  =  0pt/0o  <  1.  As  such,  these  results  should  greatly  aid 
parametric  system  modeling  of  (proposed)  LADARs  that  are  employed  for  both  target 
ranging  and  analytical  cost  engineering  modeling  of  corresponding  space  defense  systems. 
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Table  1 .  Definitions  of  Some  Key  Parameters  and  Reference  to  Their  Appearance  in  the  Body 
of  this  Report.  The  quantities  q,  t,  0pt,  and  oj  are  the  charge  on  the  electron,  return-signal 
pulse-width,  angular  boresight  evror,  and  single-axis  angular  jitter  standard  deviation, 
respectively. 


EQUATION  NO./ 

PARAMETER 

VALUE 

COMMENTS 

Ks 

Mean  signal  photoelectron 
count  per  pulse 

2 

's 

Mean  signal  current  per  pulse 

qKg/T 

m 

Mean  number  of  speckles 
contained  within  the  collecting 
aperture 

7 

6o 

Laser  beam  angular  radius 

4 

Non-signal  related  noise  current 

16 

variance 

Pfa 

Probability  of  false  alarm 

17 

ij 

Threshold  current 

Obtained  from  Eq.  (19) 

Go/ai 

11b 

b 

11c 

£m 

•t 

m 

‘s 

28 

Kd 

Mean  photoelectron  count, 
averaged  over  the  pointing 
statistics 

14 

Pd 

Detection  probability  per  pulse, 
averaged  over  the  pointing 
statistics 

12 
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Table  2.  Probability  of  Detection  per  Pulse,  Averaged  over  the  Pointing  Statistics  in  the 
High  Photoelectron  Count  Regime.  The  quantity  a  is  given  by  Eq.  (41). 


m 

Pd 

COMMENTS 

«  1 

(unresolved 

targets) 

exp  [-  ij/aisl 

Valid  for  all  values  of  noise. 

1  <  m  <  5-10 
(partially  resolved 
targets) 

T(m,€„la) 

T(m) 

Valid  lor  (SNR)q  >1.  For  integer 
m,  use  Eq.  (30)  with  replaced 
by  £fY,/a-  For  noninteger  values, 
smoothly  interpolate  between 
successive  integers 

m  >  10 

(resolved  targets) 

1 

2| 

1-Herff 

1  )\ 

Valid  for  all  values  of  Use  the 

rational  approximation,  given  in 

Sec.  IIIB,  to  the  error  function  in 
spreadsheet  programs.  For  a  given 

Op  and  Pfg,  the  threshold  current  is 
obtained  from  Eqs.  (18)  and  (19). 

where 

f  X  m 

2 

9 
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Appendix  A 
Pointing-error  Statistics 

In  this  Appendix,  we  outline  the  derivation  of  Eq.  (11).  We  assume  circularly  symmetric 
Gaussian  angular  jitter  statistics  with  a  fixed  angular  offset  bias  for  the  LADAR's  pointing 
system.  Let  0  denote  the  angle  with  respect  to  the  direction  of  the  target  at  which  the  beam 
points  at  any  instant  of  time.  The  probability  density  of  0  being  in  the  range  0  to  0  -I-  d0  is 
given  by  the  Rice-Nakagami  distribution.^  In  the  notation  used  here  this  distribution  is  given 
by 


O’. 


20-? 


.  ) 


(A-l) 


v/here  Iq  is  the  modified  Bessel  function  of  order  zero,  Oj  is  the  single-axis  angular  jitter 
standard  deviation,  and  0pt  is  the  corresponding  offset  bias.  The  corresponding  probability 
density  of  irradiance  on  the  target,  Ppt(I),  can  be  obtained  from  Ppt(6)  through  the  relation 


Pp.{I)  = 


dl 


Pp.(^) 


(A-2) 


For  a  given  0,  the  irradiance  at  the  target  location  is  given 

I(«)  =  I.  exp[- «=/2e;]  , 


from  which  it  follows  that 


A-l 


(A-3) 


Differentiating  Eq.  (A-3)  with  respect  to  I  and  substituting  this  result  together  with  Eq.  (A-1) 
into  Eq.  (A-2)  yields,  after  algebraic  simplification, 


a'e 


■■ni 


il/iy-'j,[a\l2en(lll)).  0</<, 


(A-4) 


=  0,  Otherwise  , 


where  Jq  is  the  Bessel  function  of  the  first  kind  of  order  zero,  a  =  0o/Oj,  and  b  =  Gpt/eo-  By 
noting  that  the  signal  photoelectron  count  (and  corresponding  signal  current)  is  directly 
proportional  to  the  target  illumination,  Eq.  (11)  follows  directly  from  Eq.  (A-4). 

Although  in  general  no  closed  form  integral  exists  for  Eq.  (12),  we  can,  however,  obtain  an 
analytic  result  for  the  special  case  of  no  boresight  pointing  error  (i.e.,  b  =  0).  With  Pj  given 
by  Eq.  (9),  it  can  be  shown  that  the  probability  of  detection  per  pulse,  averaged  over  the  jitter 
distribution,  is  given  by 


Pd  =1- 2^1  [a^m.  l  +  a’.-K^/m]  ,  (A-5) 

where  2F1  (a,  b,  c,  :z)  is  the  hypergeometric  function  (discussed  in  Chapter  15  of 
Reference  6).  Numerical  results,  based  on  Eq.  (A-5),  are  shown  in  Figure  4  and  7.  In  the 
limit  m  ->  00  (i.e.,  Poisson  statistics)  Eq.  (A-5)  yields  Pj  =  1  -  ,F,  (a‘ ,  1  +  a^  -  kJ,  where 

iFj  (a,b;z)  is  the  confluent  hypergeometric  function  (see  Sec,  9.2  of  Ref.  5). 
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Appendix  B 

The  Negative^binomial  Distribution: 

Some  Limiting  Cases 

First  consider  the  case  of  large  Kg  and  m.  Substituting  Stirling's  approximation  (i.e., 
r (z  +1)  =  V2^z*'^^e'^)  for  the  gamma  functions  appearing  in  Eq.  (6),  we  obtain 


PAK)  = 


_L_^/w 


(B-1) 


where 


f(K)  =  (K  +  m  -  j)  t^n(K  +  m)  -  (K  +  i)fn(K)  -{m-  \)Cn{m) 


-K^n 


+  m^ 

K.  J 


m^nf 


+m 
V  m 


(B-2) 


For  large  Kg,  we  expect  that  P(K)  is  peaked  near  K  ~  Kg.  Indeed,  differentiating  f(K)  with 


respect  to  K  yields  f  (K)  = 
neglected. 


^  1  +  m/K  ^ 


1  +  m/K 


,  where  terms  of  the  order  1/K,  1/m,  and  higher  are 


s  y 


From  the  condition  f  '(K)  =  0  we  obtain  that  the  maximum  value  of  f  is  achieved  at 


K  =  K. 


(B-3) 


B-1 


Thus,  expanding  f(K)  in  a  Taylor’s  series  about  Kg  yields  to  terms  of  second  order 


f(K)  =  /(AT,)  +(*•  -  K,)r(K,)*  UK  -  K,)'r(K,) 


-  \£na\ 


2ctI 


(B-4) 


where 


K: 


m 


(B-5) 


Hence,  from  Eq.  (B-4)  and  Eq.  (B-1),  we  obtain 


P,(K)  =  -pi~exp 


(K-K,)^  I 

J 


(B-tS) 


Thus,  for  large  Kg  and  m,  the  negative-binomial  distribution  is  Gaussian,  with  the  same  mean  and 
variance  of  the  underlying  distribution.  The  corresponding  probability  distribution  of  signal 
current,  Eq.  (20),  follows  directly  from  Eq.  (B-6)  since  i  =  qK/x.  In  particular,  we  have 
a:  ^  Var[i]-(q/r)'Var(K]=  (q/r)'  (k.  +K:  /m)  =  qi,/r+  i'/m,  where  the  mean  signal 

current  ig  =  qKg/x. 


An  elementary  argument  involving  the  Central  Limit  theorem  can  be  employed  to  derive 
Eq.  (B-6)  directly.  Each  of  the  m  independent  correlation  cells  or  speckles  contained  with  the 


collecting  aperture  leads  to  a  photocount  distribution  that  obeys  Bose-Einstein  statistics  with  a 
mean  count  kj  =  Kg/m,  and  a  variance  =  kg  +  kg  =  Kg/m  +  (Kg/m)^ .  The  total 

photocoupt  from  the  entire  collecting  aperture  is  just  the  sum  of  the  (non-zero)  contributions 
t^rom  each  of  the  m  independent  correlation  cells.  The  photocount,  k,  due  to  a  single 
correlation  cell  is  of  the  order  K/m.  Hence  for  k  to  be  non-zero,  it  follows  that  K  must  be 
large  since  m  >>  1.  Thus,  for  both  large  m  and  K  we  have,  by  the  Central  Limit  theorem,  the 
resulting  probability  distribution  approaches  the  Gaussian  distribution  with  a  mean  equal  to 
mkg  =  Kg  and  a  variance  equal  to  mo^  =  Kg  +  Kg /m ,  in  agreement  with  Eq.  (B-6). 

Next  consider  the  negative-binomial  distribution  in  the  limit  m  — >  1 .  In  this  case,  we  obtain 
from  Eq.  (6) 


P,(K)  = 


r  K. 


m  =  1 


(B-7) 


This  is  known  as  the  Bose-Einstein  (or  geometric)  distribution.  For  large  photoelectron 
counts,  we  have  (Kg/Kg  -I-  1)*^  =■  exp(-  K/KJ,  and,  hence,  Eq.  (B-7)  becomes 


P,(K)^ 


exp(-  K /  AT,) 


K, 


(B-8) 


which  is  known  as  the  Rayleigh  (or  negative-exponential)  distribution. 

On  the  other  hand,  consider  the  limit  of  the  negative-binomial  distribution  as  m  x.  From 
Eq.  (6),  we  have 


B-3 


P,(K)  = 


1  lim 

K!  m  ->  00 


r(K  +  m) 


r(m) 


1  + 


m 


TT  l-m 

1  +  ^ 
m 


1  K^K  .  _K 

—  —  e 

K!  Vm 


K„-K. 


K.^^e 


K! 


(B-9) 


which  is  the  Poisson  distribution. 

Finally,  we  show  that  for  large  K,  the  Poisson  distribution  becomes  a  Gaussian  distribution 
with  a  mean  equal  to  and  a  variance  equal  to  the  mean.  Of  course,  this  result  follows 
directly  from  Eq.  (B-6).  However,  since  the  Poisson  distribution  is  interesting  in  its  own  right 
and  plays  such  an  important  role  in  probability  and  statistics,  we  give  an  independent 
d.mvation,  that  starts  directly  from  the  Poisson  distribution  itself.  For  large  K,  we  may  use 
Stirling's  approximation  for  K!  and  express  the  Poisson  distribution  as 


P(A')  = 


(B-10) 


where 


g(K)  =  K  +  K/:nK, -(K+^)fnK-K,  (B  11) 

This  function  has  a  stationary  point,  for  large  K,  at  K  ==  K^.  Expanding  g(K)  in  a  second  order 
Taylor  series  about  K  =  yields 


B-4 


(B-12) 


Hence,  substituting  Eq.  (B-12)  into  Eq.  (B-10)  yields 


TO  = 


2K, 


(B-13) 


which  is  the  Gaussian  distribution  with  the  variance  equal  to  the  mean. 


B-5 


Appendix  C 

Signal-noise-limited  Detection  Statistics 


In  this  appendix,  we  consider  signal-noise-limited  detection  statistics  where  the  mean  number  of 
signal  photoelectrons  per  pulse  »  1 ,  and  the  mean  number  of  non-signal-related  noise  photo- 
electrons  is  negligible.  The  negative-binomial  distribution  given  by  Eq.  (6)  is 


P.(K)  = 


r(m)r(K)  I,  kJ  I  mj 


(C-l) 


For  Kg  »  1,  PgCK)  is  nonzero  for  K  large.  In  chis  case,  we  have 

raC  +  m- l)  =  Km*l  r(K)  . 


and 


-K 

mK' 

1  +  — 

1  kJ 

=  exp 

.  K. 

Thus 


1  1 

m 

r(m)l 

K” 


exp 


(C-2) 


Next  we  show  that  for  large  Kg  and  m  the  detection  probability,  Pj,  given  by  Eq.  (27)  is  identical 
to  that  obtained  from  Eq.  (24)  in  limit  n  — »  0.  From  Eq.  (27),  we  have 


C-l 


in£ 

For  large  m,  the  integrand  in  Eq.  (C-3)  is  peaked  about  t  =  m  -  1.  We  now  approximate  the 
integrand  in  a  particular  way,  by  writing  it  in  an  exponential  form,  exp  [f(t)],  and  using  a  Taylor 
series  approximation  for  f(t)  near  its  maximum.  (A  Taylor's  expansion  of  the  integrand  itself, 
for  example,  would  not  be  useful  if  only  a  few  terms  are  kept.)  The  integrand  of  Eq.  (C-3)  is 

e^“)  =  t'"'‘e“' =exp[(m  -  l)£nt  - 1]  , 

so 

f(t)=(m-l)^nt-t 

f'(t)  =  Hl — f'=0  for  t  =  .m-l 
t 

(m-1) 
t^  ’ 

then,  expanding  about  t  =  m  -  1  yields 


C-2 


P  (^-0'"~*exp[-(m-l)]7 

"  r(m)  J 


exp 


mt 


dt 


^n/2(m-\)(m  -  if  exp[-(m- 1)] 
_  _ 


l+erf 


Vm-l(l-e) 


(C-4) 


For  large  ra,  r(m)  ==  ^^ex  Substituting  this  into  Eq.  (C-4)  yields 


r 


I + erf  I 


Vm(l-e) 

~w~ 


(C-5) 


which  is  identical  lo  Eq.  (24)  for  Osn  -Cfs  =  /™-  Thus,  in  the  limit  of  large  m,  Eq.  (24)  and 

Eq.  (27)  are  identical  for  signal-noise-limited  detection.  In  Figure  C-1,  we  have  compared  these 
two  expressions  for  m  =  5  and  10.  For  engineering  purposes,  use  of  Eq.  (24)  for  m  >  5  should 
be  adequate.  Incidentally,  it  follows  directly  from  Eq.  (C-3)  and  (C-5)  that 


lim  r(m,me) 
m— r(/n) 


1. 

0<e  <  r 

e  =1 

0, 

E  >  1 

(C-6) 


This  interesting  result,  to  the  best  of  my  knowledge,  has  not  yet  appeared  explicitly  in  the 
literature. 


C-3 


m  =5 10 


Figure  C-1.  Signal-noise-limited  detection:  a  comparison  of  tlte  detection  probability  as  a 

function  of  the  signal-to-threshold  ratio.  The  dashed  curve  is  based  on  Eq.  (24)  and 
the  solid  curve  is  based  on  Eq.  (27). 
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Appendix  D 

Probability  Distribution  of  a  Signal  Photoelectron  Count  Obeying  Bose- 
Einstein  Statistics  Immersed  in  Poisson  Noise 


For  m  =  1,  Eq.  (25)  becomes 


l  +  K„ 


1  + 


(D-1) 


where 


3=n(l  +  K,-') 


(D-2) 


The  sum  appearing  in  Eq.  (D-1)  exists  in  closed  form.^  The  result  is 


y  r(ir-fl.p) 

^  r(/i:  +  i) 


(D-3) 


where  r(a^)  is  the  complementary  incomplete  gamma  function.  Substituting  Eq.  (D  3)  into 
Eq.  (D-1)  yields 


l  +  K  AT!  •  . 


(D.4) 


which  is  valid  for  arbitrary  value  of  the  mean  photoelectron  count,  Kj,  and  non-signal-related 
mean  photoelectron  count,  n.  It  can  be  shovra  for  arbitrary'  Kg  and  n  that  E[K]  =  Kg  +  n  and  Var 
[K]  =  Ks+K2  +  n. 
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In  Figures  D-1  through  D-4,  we  plot  as  a  function  of  K  for  a  representative  range  of  values 
of  Kj  and  n.  Except  for  large  photocounts,  no  further  simplication  of  Eq.  (D-4)  has  been 
obtained.  For  Kj  »  1,  Eq.  (D-4)  can  be  written  as 


P„(K) 


e,p  [-(K-n)/K.]^(^  „) 

K. 


(D-5) 


where 


G(K,n) 


r(K-i-l,n) 

r(K  + 1) 


(D-6) 


Following  the  procedure  outlines  in  Appendix  B,  it  can  be  shown  for  large  K  that 


0,  K  <  n 
M,  K  =  n 
1,  K>n 


(D-7) 


That  is,  G  behaves  as  a  unit  step-function  located  at  K  -  n.  These  features  are  illustrated  in 
Figure  D-5  for  K  ~  1000.  Hence, 


(D-8) 


D-2 


K 

Figure  D-1.  Probability  density  of  a  signal  obeying  Bose-Einstein  statistics  immersed  in  Poisson 
noise,  where  the  mean  signal  count  =  1  for  various  values  of  the  noise  count.  For 
convenience  of  presentation,  a  continuous  curve  is  drawn  between  integer  count 
values. 
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Figure  D-2.  Probability  density  of  a  signal  obeying  Bose-Einstein  statistics  immersed  in  Poisson 
noise,  where  the  mean  signal  count  =  10  for  various  values  of  the  noise  count. 
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Figure  D-3.  Probability  density  of  a  signal  obeying  Bose-Einstein  statistics  immersed  in  Poisson 
noise,  where  the  mean  signal  count  =  100  for  various  values  of  the  noise  count. 
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Figure  D-4.  Probability  density  of  a  Signal  obeying  Bose-Einstein  statistics  immersed  in  Poisson 
noise,  where  the  mean  signal  count  =  1000  for  various  values  of  the  noise  count. 
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Figure  D-5.  Quantity  G(K,  n)  =  r(K  +  1, 
values  of  n. 
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where  U  is  the  unit  step-function.  Thus,  in  this  Unit,  the  probability  of  signal  plus  noise  photo¬ 
electron  counts  is  a  shifted  Bose-Einstein  distribution,  the  shift  being  equal  to  the  mean  noise 
photocount.  The  corresponding  probabiUty  of  detection  is  given  by 


f^T 

=  ]  P.WdK 

ii+AATj- 


(D-9) 


where  AKj  is  the  threshold  co  mt  relative  to  the  noise  floor  count,  n.  For  large  counts,  where 

one  oper^aonally  measures  a  current  rather  than  discrete  electron  counts,  the  probabiUty  of 
detection  is  given  by  exp(-  i  i  /ij),  where  ij  is  the  threshold  current  relative  to  the  mean  non¬ 
signal  related  noise  cur 'em,  lu,  and  ij  is  the  corresponding  mean  signal  current.  Unfortunately, 
other  than  for  m  -  1,  no  analytic  closed-form  expressions  have  been  obiair  ed  for  Eq.  (25). 

Finally,  consider  he  limit  of  the  probability  of  emission  of  K  signal  plus  loise  photoelectrons, 
where  the  mean  '  ignal  photocount  is  much  less  than  the  corresponding  noise  photocount  (i.e., 

Kj «  n).  In  this  case  we  may  use  the  asymptotic  expansion  of  the  complementary  incomplete 
gamma  function^, 

r(a,z)»e-^z‘’-'|^l  +  ^  +  ...j  ,  (D-10) 


0-6 


in  Eq.  (D-4)  to  obtain 


_  (1+C‘r  ,,K.  n^(\+K:'f(l  +  KKJn) 

"  1+a:.  a:! 


^  1  V'"ii*^(i+Aa:,/w) 

\^J 


AT! 


^  exp[-(/i+A:,)](/i+A:,/ 
AT! 


to  first  order  in  Kj.  Thus,  for  Kg  «  n  and  arbitrary  values  of  n,  the  statistics  of  the  total 
emission  process  are  Poisson. 
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TECHNOLOGY  OPERATIONS 


Tne  Aerospace  Corporation  functions  as  an  "architect-engineer"  for  national  security 
programs,  specializing  in  advanced  military  space  systems.  The  Corporation's  Technology 
Operations  supports  the  effective  and  timely  development  and  operation  of  national  security 
systems  through  scientific  research  and  the  application  of  advanced  technology.  Vital  to  the 
success  of  the  Corporation  is  the  technical  staffs  wide-ranging  experti.'e  and  its  ability  to  stay 
abreast  of  new  technological  developments  and  program  support  issues  associated  with  rapidly 
evolving  .space  systems.  Contributing  capabilities  are  provided  by  the.se  individual  Technology 
Centers: 

Electronics  Technology  Center:  Microelectronics,  solid-state  device  physics, 

VLSI  reliability,  compound  semiconductors,  radiation  hardening,  data  storage 
technologies,  inf'ared  detector  devices  and  testing;  electro-optics,  quantum  elecnonics, 
solid-stdte  lasers,  optical  propagation  and  communications;  cw  .and  pulsed  chemical 
laser  development,  optici  resonators,  beam  control,  atmospheric  propagation,  and 
laser  effects  and  countermeasures;  atomic  frequency  standards,  applied  laser 
spectroscopy,  laser  chemistry,  laser  optoelectronics,  pha^  conjugation  and  coherent 
imaging,  solar  cell  physics,  battery  electrochemistry,  battery  testing  and  evaluation. 

Mechanics  and  Materials  Technology  Center:  Evaluation  and  characterization  of 
new  materials:  metals,  alloys,  ceramics,  polymers  and  their  composites,  and  new 
forms  of  carbon;  development  and  analysis  of  thin  films  and  deposition  techniques; 
riondesUTictive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  coirosion;  development  and  evaluation  of  hardened  components; 
analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures;  launch 
vehicle  and  reentiy  fluid  mechanics,  heat  u-ansfer  and  flight  dynamics;  chemical  and 
electric  propulsion;  spacecraft  structural  mechanics,  spacecraft  survivability  and 
vulnerability  assessment:  contamination,  thermal  and  structural  control;  high 
temperature  thermomechanics,  gas  kinetics  and  radiation;  lubrication  and  siiriace 
phenomena. 


Space  and  Environment  Technology  Center:  Magnctospheric,  auroral  and 
cosmic  ray  physic,s,  wave-particle  interactions,  magnctospheric  plasma  waves; 
atmospheric  and  ionospheric  physics,  density  and  composition  of  the  upper 
atmosphere,  remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared 
astronomy,  infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere;  effects 
of  electromagnetic  and  paniculate  radiations  on  space  systems;  space  instrumentation; 
propellant  chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
atmo,spheric  chemical  reactions,  atmospheric  optics,  light  .scattering,  state-specific 
chemical  reactions  and  radiative  signatures  of  missile  plumes,  itnd  sensor  out-of-field- 
of-view  rejection. 


